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Visualization of chiral condensate at finite temperature on the lattice 
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We perform an analysis of the topological and chiral vacuum structure of four-dimensional QCD on the lattice at 
finite temperature. From correlation functions we show the existence of local correlations between the topological 
charge density and the quark condensate on gauge average. We comment on sizes of clusters of nontrivial chiral 
condensate and of instantons in full QCD. By analysis of individual gauge configurations, we demonstrate that 
at the places in Euclidian space-time, where instantons are present, amplified production of quark condensate 
occurs. 



It is believed that the instantons as carriers of 
the topological charge might play a crucial role 
in understanding the confinement mechanism of 
four-dimensional QCD, if one assumes that they 
form a so-called instanton liquid Instan- 
tons have topological charge Q being related to 
the zero eigenvalues of the fermionic matrix of 
a gauge field configuration via the Atiyah-Singer 
index theorem ||^ . Recently, it was demonstrated 
that monopole currents which constitute a dif- 
ferent topological excitation of compact SU(3) 
gauge theory, appear preferably in the regions of 
non- vanishing topological charge density . It 
has been conjectured that both instantons and 
monopoles are related to chiral symmetry break- 
ing ||l|,||j6j. In this contribution we further sup- 
port this idea by the following results of a direct 
investigation of the local correlations of the quark 
condensate and the topological charge density on 
realistic gauge field configurations. 

For the implementation of the topologi- 
cal charge on a Euclidian lattice we re- 
strict ourselves to the so-called field theo- 
retic definitions which approximate the topolog- 
ical charge density in the continuum, q{x) = 
^e'^'^P'^ Tr (^Ff,^{x)Fp^{x)j . We used the pla- 
quette and the hypercube prescription. To get rid 
of quantum fluctuations and renormalization con- 
stants, we employed the Cabbibo-Marinari cool- 
ing method. Mathematically and numerically the 
local chiral condensate 'iljip{x) is a diagonal ele- 



ment of the inverse of the fermionic matrix of the 
QCD action. We compute correlation functions 
between two observables Oi{x) and 02(2/) 
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and normalize them to the smallest lattice separa- 
tion dmin, c{y-x) = g{y-x)/g (drain)- Since topo- 
logical objects with opposite sign are equally dis- 
tributed, we correlate the quark-antiquark den- 
sity with the square of the topological charge 
density. Our simulations were performed for 
full SU(3) QCD on an 8^ x 4 lattice with pe- 
riodic boundary conditions. Applying a stan- 
dard Metropolis algorithm has the advantage that 
tunneling between sectors of different topologi- 
cal charges occurs at reasonable rates. Dynami- 
cal quarks in Kogut-Susskind discretization with 
Uf = 3 flavors of degenerate mass m = 0.1 were 
taken into account using the pseudofermionic 
method. We performed runs in the confinement 
phase at /3 = 5.2. Measurements were taken on 
1000 configurations separated by 50 sweeps. 

Figure |l| shows results for the correlation 
functions of Eq. (0) with Oi = -tp^p^x) and 
^2 = I'^iy) (plaquette definition). The -ipipq^- 
correlation function exhibits an extension of more 
than two lattice spacings, indicating nontrivial 
correlations. To gain information about the cor- 
relation lengths exponential fits to the tails of the 
correlation functions were performed. They show 
that an increasing number of cooling steps yields 
shorter correlation lengths. The corresponding 
screening masses are C, = 0.59 and C — 1.56 in 
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Figure 1. Correlation function of the quark- 
antiquark density and the topological charge den- 
sity for and 11 cooling steps. The correlations 
extend over two lattice spacings and indicate local 
coexistence of the quark condensate and topolog- 
ical objects. 

inverse lattice units for and 11 cooling steps, 
respectively. They have to be interpreted as ef- 
fective masses and reflect the effective gluon ex- 
change and the vacuum structure of QCD. It is 
assumed that the size p of a t'Hooft instanton 
qp{x) ~ yo'^ (x^ -I- yO^)~'* centered around the ori- 
gin enters also into the associated distribution of 
the chiral condensate ^iljp{x) ~ (x^ -I- p^)~^ 
0. To estimate p we fitted a convolution of the 
functional form /(a;) = J ijjijjp{t)qp{x~t) dt to the 
data points. This was evaluated after 11 coohng 
steps where the configurations are reasonably di- 
lute. Our fit yields p(V'V"Z^) — 1-8 in lattice spac- 
ings. To check consistency we extracted from the 
gg-correlation a value of p{qq) = 2.05 which is in 
good agreement. 

One may choose to directly visualize densities 
of the quark condensate and topological quanti- 
ties on individual gauge fields rather than per- 
forming gauge averages. We persue this in the 
following to get first insight into the local inter- 
play of topology with chirally nontrivial regions. 
In a series of papers we found that at the local 
regions of clusters of topological charge density, 
which are identified with instantons, there are 
monopole trajectories looping around in almost 
all cases for both pure SU(2) and SU(3) gauge 
theories In Fig. ^ three typical topologically 
nontrivial configurations from SU(3) theory with 



dynamical quarks on the 8'^ x 4 lattice in the con- 
finement phase are shown for fixed time slices. 
We display the topological charge density (hyper- 
cube definition) by dark dots if the absolute value 
\q{x)\ > 0.003. The quark-antiquark density is 
indicated by light dots whenever a threshold for 
i>il}{x) > 0.066 is exceeded. Monopole currents 
are defined in the maximum Abelian projection 
and only one type is shown by lines. By analyz- 
ing dozens of gluon and quark field configurations 
we found the following results. The topological 
charge is hidden in quantum fluctuations and be- 
comes visible by cooling of the gauge fields. For 
cooling steps no structure can be seen in g(a;), 
ip'ip{x) or the monopole currents, which does not 
mean the absence of correlations between them. 
After about 5 cooling steps clusters of nonzero 
topological charge density and quark condensate 
are resolved. These particular configurations con- 
tain a single instanton, a pair of antiinstantons 
and an instanton-antiinstanton pair. For more 
than 10 coohng steps both topological charge and 
chiral condensate begin to die out and eventually 
vanish. Combining the above finding of Fig. 1 
showing that the correlation functions between 
ipipi^) s-nd q^{y) are not very sensitive to cool- 
ing together with the cooling history of the 3D 
images in Fig. 2, we conclude that instantons go 
hand in hand with clusters of ■tptpix) ^ also in 
the uncooled QCD vacuum 

Figure 2 has demonstrated that the quark- 
antiquark density attains its maximum values at 
the same positions where the extreme values of 
the topological charge density are situated. This 
behavior is further substantiated in Fig. ^ where 
the ?/5-(/'(a;)-values are plotted against q{x) for all 
points X. The nine selected configurations at 10 
cooling steps have different topological content. 
At first sight a linear relationship between the 
absolute value of the topological charge density 
and the virtual quark density is suggested. 

We comment on the behavior of the topological 
structure and the quark condensate when crossing 
the phase transition. The normalized correlation 
functions do not change qualitatively The 
strength of the correlations becomes, however, 2- 
3 orders of magnitude smaller. This means that 
the topological activity becomes much weaker in 
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Figure 2. Cooling history for time slices of three gauge field configurations of SU(3) theory with 
dynamical quarks. The columns represent a single instanton, a pair of antiinstantons and an instanton- 
antiinstanton pair as the configuration gets cooled. The dark dots represent the positive and negative 
topological charge density; the light dots the density of the quark condensate, ft turns out that the quark 
condensate takes a non- vanishing value at the positions of instantons. 
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Figure 3. Scatter plots of the local chiral condensate ^tp{x) and the topological charge density q{x) for 
nine configurations after 10 cooling steps. A linear relationship is suggested. 
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the deconfinement phase as expected. Most of the 
configurations have trivial net topological charge. 
This does not exclude the existence of instanton- 
antiinstanton pairs which become more difficult 
to be resolved at the smaller physical volume at 
f3 = 5.4 which we also considered. In less than 
one percent of the gauge field configurations in- 
stantons could be identified. The configurations 
we scanned did not show pronounced pictures of 
instantons and quark condensate. The quark- 
antiquark density is considerably lower and does 
not have the tendency to cluster anymore, ipip^x) 
becomes uniformly distributed over the lattice 
sites. The maximum values of il}ip{x) in configura- 
tions in the deconfinement after 10 cooling steps 
are one order of magnitude smaller. 

In summary, our calculations of correlation 
functions between topological charge and the 
quark condensate yield an extension of about two 
lattice spacings. The correlations suggest that 
the local chiral condensate takes a non-vanishing 
value predominantly in the regions of instantons 
and monopole loops. It was well known before 
that the chiral condensate is related to the topo- 
logical charge and topological susceptibility. The 
visualization exhibited that the distribution of 
the "chiral condensate" concentrates around ar- 
eas with enhanced topological activity (instan- 
tons, monopoles). We demonstrated that exactly 
at those places in Euclidian space-time, where 
tunneling between the vacua occurs, amplified 
production of quark condensate takes place. It 
must be emphasized that this represents the sit- 
uation on a finite lattice with finite quark mass 
without the extrapolation to the thermodynamic 
and chiral limit. We found for full SU(3) QCD 
with dynamical quarks that the clusters of non- 
vanishing quark condensate have a size of about 
0.4 fm, which corresponds to the instanton sizes 
observed in the same configurations. Visualiza- 
tion of quark and gluon fields might be especially 
useful to decide if the instanton-liquid model is 
realized in nature and if instanton-antiinstanton 
pairs appear in the deconfined phase. It might 
also help to clarify the question of the existence of 
a disoriented chiral condensate with consequences 
for heavy-ion experiments. 
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DISCUSSIONS 

Heinz J. Rothe, Univ. Heidelberg 

1. Your results suggest that the dynamical 
mechanism for the deconfinement phase 
transition and chiral symmetry restauration 
is the same. Is this also your opinion? 

2. How does the cooling procedure affect the 
number of monopoles you see? 

Harald Markum 

1. It was shown by lattice simulations that the 
temperature and chiral phase transition co- 
incide. My opinion is that there is a com- 
mon driving force, maybe the gluon field it- 
self 

2. Cooling smooths quantum fluctuations and 
uncovers topological charges. One of the 
caveats is that monopoles are gradually lost. 
This is a reason why one tries to develop al- 
ternative methods. 



